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摘 E: ARM Magnus FER MIRE Schrödinger Zz. Schrödinger HFAA YY FIERE, 
用 适当 差分 格式 对 其 进行 模 平 方 守恒 空间 离散 ， 转 化 成 模 平 方 守恒 的 常 微分 方程 组 ， 再 
用 Magnus 方法 求解 。 数 值 结 果 表 明 Magnus 方法 能 保 非 线性 Schrödinger 方程 模 守 恒 量 的 优越 
性 和 好 的 稳定 性 。Magnus 方法 可 应 用 到 其 它 模 守恒 的 微分 方程 。 
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1 引言 


近年 来 对 数学 和 物理 中 的 微分 方程 数值 模拟 时 ， 不 仅 要 能 正确 地 模拟 微分 方程 的 变化 行 
为 ， 而 且 要 使 数值 解 能 尽量 地 保持 微分 方程 的 某 些 固有 特性 ， 如 非 线性 Schrödinger 方程 


ide 十 Wzz 十 alla = 0, (1) 


初始 条 件 V(z,0) = polz) CER. a> 0 是 常数 。 方 程 (1 上) 是 孤立 子 理论 中 最 重要 的 可 积 模 
H, 广泛 应 用 在 物 埋 中 的 许多 领域 ， 包 括 非 线性 光学 和 等 离子 体 物理 。 方 程 (1) 具有 守恒 特性 


gal “Wed ae = Qo, (2) 


emf (3f + fur ]ee= a 


其 中 Qo 和 Eo 是 常数 。Sanz-Serna 提出 了 蛙 跳 格 式 和 改进 的 Crank-Nilcolson 格式 求解 Schréd- 
inger 方 程 ， 讨 论 了 这 些 格式 是 否 保 方程 的 守恒 特性 趾 。 在 文献 [2,3] 中 ， 提 出 了 新 的 模 守 恒 差 
分 格式 求解 Schr6dinger 方程。 然而 常见 的 差分 格式 一 般 都 是 隐 式 的 ， 精 度 上 低 阶 。 

近年 来 在 计算 数学 领域 内 产生 了 一 个 新 的 分 支 一 保 结构 算法 。 该 方法 不 仅 具 有 经 典 算法 
的 高 阶 精度 ， 稳 定性 ， 健 壮 性 ， 长 时 间 精 确 模拟 微分 方程 的 变化 行为 ， 并 能 保 微分 方程 的 一 
些 内 部 守恒 特性 ， 共 中 包括 冯 康 和 其 保 结构 算法 研究 组 成 员 发 展 的 辛 算法 ， 在 准 算 法 基础 上 
提出 的 多 学 算法 ， 和 挪威 的 Munthe-Kaas， 英 国 的 Iserles 等 人 发 展 了 保 一 般 李 群 方程 不 变性 
的 李 群 方法 则 ,45 。 李 和 群 方法 包括 常见 的 RKMK 方法 ，Magnus 方 法 ，Fer 方法 等 。 在 这 里 ， 利 
用 李 和 群 方法 中 的 Magnus 方 法 构造 了 模 平方 守恒 格式 ， 用 Runge-Kutta 方 法 和 Magnus 方法 求 
解 Schrödinger JFT, RANEE BIR Schrödinger 方程 的 守恒 特性 。 
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2 非 线 性 Schr6dinger 方程 的 离散 
Wy =ptig, JERE Schrodinger 77 ## (1) 等 价 于 
pt + qzz + a(p” + q7)q = 0, (4) 


dt — Prz — alp? + q’)p = 0. (5) 


/ "i 


的 守恒 特性 。 数 值 离散 时 ， 要 使 离散 格式 尽量 保持 方程 变量 的 离散 守恒 ， 能 反映 方程 (1) 的 守 
恒 特 性 。 
方程 (4)-(5) 变 为 


方程 (1) 具有 


其 中 


( 
是 一 个 微分 算 子 ，c = a(p? + 92) 和 了 是 单位 算 子 。 对 方程 (6) 进行 适当 空间 离散 。 设 对 eB 取 
二 阶 中 心 差分 ， 方 程 (6) 离散 成 


“Ti 1 0 0 
dz 1 0 一 4 1 —T2 1 0 
= as 和 ， (7) 
A 0 
0 0 0 TN 


其 中 

Tm =2— alp + Ar, k=l2N, z= (phpa PN Gas" dN) 
i dz = A(z)z, A(z) RR EMM. HFE (6) MEA Ol = X% | Raa 
3 Magnus 方法 

设 流 形 M 上 的 常 微分 方程 
Y =A(t)Y¥, t>0, Y(O)EG, (8) 

HPA: R > g fll A(t) ERM, ACGMY EG, (GRATE, g 为 相应 的 李 代数 )。 
方程 (8) 的 解 为 Y(t) = exp (O(t)) Yo, t > 0， 其 中 


日 =dexp5l4=》， Zt adh A, t>0, ©(0)=0, (9) 
k=0 ` 
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{Br}kez+ Æ Bernoulli RX. FA Picard E(t FR (9), 2 (t) = O 


[m-+1] ' =i S Bef” jk 
© o= | dexpain AO = °F J) ad@tm(e) A(E)dE, m=0,1,--- (10) 
tn O(t) = fy A(&1)d&i 和 
aBa = [anes 3/ pi 4(6ajdEz, A(E) | dE 
Lf agoae[ { Aeae aa + a1) 
由 Picard 定 理 可 知 6( = lim OM (t) 4) Hea LBL A — Ah) SB BC HE, ERE O(t) fé A HE 
阵 4 ARDIA FRET IE, BI Magnus 展开 
= >》 H(t), 
k=0 
Hi 足 一 个 包含 十 1 阶 重 积分 的 线性 组 合 ， 如 
Ho(t) = 人 A(&1)déi, me =—5 [| if A(é2)dé2, 4(6)] da. 
对 Magnus 展开 中 的 积分 形式 进行 数值 积分 ， 取 
t t pi 
= | A9 ny - be | (ae). a@n)ae, 
s0= {ff aeaea] 2 S LP (ve Ae) ey) 
每 一 个 Magnus 展开 项 表示 为 了 ) = fs L(A(Er), A(é2), + , A(Es))d€, 工 是 一 个 多 项 式 ， Skt 


求 积 区 间 内 离散 点 的 集合 
S= {été€R’:é € [0,h], & € [0,m,], l= 2,3,... ,s), 


m € {1,2,3,---,l-1}, 1=2,3,---,8 


按 如 下 方式 离散 I(h)。 Hv TAD Ka €1,€2,°°' ,Cv €E [0, lj, $ Ak = hA(cxh), k= 1,2,- U 


得 积分 
K(h) = > be L(Ak,; Akas*** » Ake); 
kecy 


Cs 是 集合 {1,2,… ,v} 中 元 素 所 组 成 的 集合 ， 其 中 bs = fz [i lr (E)E 和 
$= {é E R° :& € [0, 1], & E [0, mul, l= 2,3,- sh 


以 及 


v 
— C; 
L(x) = Il z I j =1,2,°+-,v 


CG aE 
i=1,ižj * J 


273 


可 
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是 由 积分 点 构成 的 插值 多 项 式 


Ae) =i Dh) 
取代 A(t)， 并 数值 积分 ， 有 f? A(ti)dti ~ 3h(A1 + Ao), 


J [A(ti), A(t2)] dt2dt1 ~ -ltl 


hoe [ A(t), [A(t2), A(ta)|]dtadtadts = rè [ (È - FJa- (24 YB) ao, jn, Al], 
moa [A(t1), A(t2)], Alta) ]atsdtadty = A (= + sojin (=- YB) Aas lAr, Ae] 


= A((}— ¥8)h), Ap = A((4 + £h) H Alt, y) BAR A(t), FA Lagrangian 多项式 


x _ pi È 2 =1.2... 
A(t,y) =h dole (Z) Ae Ak = hA(ceh, Xz), k= 1,2,...,v, 


Ki(h) 一 > ay; L(Aji, 4j2， +: , Ajs) ak sof Thu Ei) dé, 


了 ECY Sk i=1 
Se = {€ER*:& E [0,ce], aelo Eml, 1=2,3,---, 8}. 
显 式 积分 ， 可 得 解 常 微分 方程 组 
dY 


ae TAGYY, Y(0)=¥%, AY) = —AT(t,y) (12) 


的 三 阶 显 式 Magnus H 
Ə = i(4 +442 + As) — | A1 ~ As, 242 + 2s J, 
Y?+1 — ey”, (13) 
这 里 
At =hA(tn,Y"), 42 = hA( tn + P esy”), 
Ag = hA(tn +h,e0^+4yYn), [A, B] = AB — BA. 


利用 格式 (13) 解 方程 (12)， 方 程 (12) 的 数值 解 具有 ||yn|| = [YO] KERE. 


4 ”数值 试验 


为 了 测试 显 式 Magnus 方法 能 否 保 离散 的 非 线性 Schr5dinger 方 程 的 模 守 恒 特 性 ， 定 义 模 平 
方 和 误差 


Err(t) = (¥?(0) + Y2 (0) +--- + Yov(0)) — (YPE) + Y(t) +--+ Yon(t)), t=nh, 
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初始 条 件 W%(z,0) = sech(x), z € [一 20,20]。 在 图 1 中 ， 取 a = 2， 图 1 的 右边 是 h = 0.03 时 
用 Runge-Kutta(RK) 方 法 解 方 程 (7) 得 到 的 模 平 方 和 误差 ， 误差 有 明显 增长 ， 图 1 的 左边 是 
在 疡 = 0.1 时 用 Magnus 方 法 解 方 程 (7) 得 到 的 模 平 方 和 误差 ， 误 差 仅 达 10- 蕊 ， 可 和 忽略。 在 
图 2 中 ， 图 2 的 右边 是 h = 0.375 时 用 RK 方法 解 方程 (7) 得 到 的 模 平 方 和 误差 .方程 的 模 平方 
和 趋向 无 穷 大 ， 受 步 长 限制 不 收敛 。 图 2 的 左边 是 hh = 0.8 时 用 Magnus 方 法 解 方程 (7) 得 到 的 
横 半 方 和 误差 ， 误 差 仍 仅 达 10-2。 在 很 大 步 长 下 ， 该 方法 都 稳定 。 从 图 1-2 知 ，Magnus 方 
法 能 保 Schradinger 方 程 的 模 平方 守恒 ， 而 相应 阶 显 式 RK NAT hetk Schrödinger V FERRE 
方 守恒 ， 同 时 Magnus 方 法 比 相 应 RK 方 法 在 稳定 性 方面 优越 。Magnus 可 应 用 到 类 似 模 守恒 方 
程 ， 如 Kdv 方 程 ， 耦 合 Schr6dinger 方程 等 。 


x10" 
6 0.014 
a 
| | 
-z 15 30 45 0 0 15 30 45 


t t 


图 1: M Magnus 方法 和 RK 方法 得 Schr6dinger F424 t € [0, 45] 的 模 平 方 和 误差 


225 | 45 75 


图 2: 在 不 同时 间 步 下 Magnus 方法 和 RK 方法 得 Schrödinger 方程 的 模 平方 和 误差 
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Explicit Modulus Conserving Scheme for the 
Nonlinear Schrödinger Equation 
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Abstract: This paper applies apply the Magnus method to solve the nonlinear Schrödinger equation. 
The nonlinear Schrödinger equation has the modulus conserving property. The nonlinear Schrödinger 
equation is discretizated in the spacial direction by the proper difference scheme, which is transformed 
into the ordinary differential equations. The ordinary differential equations are solved by the Magnus 
method. Numerical results show that the Magnus method have the advantage of the modulus conserving 
property of the nonlinear Schrödinger equation and its good stability. The Magnus method can be 
applied to other modulus conserving differential equations. 
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